Survival of short-range order in the Ising model on negatively curved surfaces 
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We examine the ordering behavior of the ferromagnetic Ising lattice model defined on a surface 
with a constant negative curvature. Small-sized ferromagnetic domains are observed to exist at 
temperatures far greater than the critical temperature, at which the inner core region of the lattice 
undergoes a mean-field phase transition. The survival of short-range order at such high temperatures 
can be attributed to strong boundary-spin contributions to the ordering mechanism, as a result 
of which boundary effects remain active even within the thermodynamic limit. Our results are 
consistent with the previous finding of disorder-free Griffiths phase that is stable at temperatures 
lower than the mean-field critical temperature. 
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I. INTRODUCTION 

The two-dimensional Ising lattice model is one of the 
simplest models of second-order phase transitions 0, Q • 
This model has long played a fundamental role in statisti- 
cal physics due its broad applicability and the availability 
of analytic solutions. Previous studies have proven that 
the critical behavior of this model is universal to a large 
extent, depending on the essential symmetries inherent 
to the system Q. 

The two-dimensional Ising lattice model is usually as- 
signed to a flat plane. In the last two decades, however, 
there has been a growing interest in the nature of the 
Ising model assigned to curved surfaces [3, @, 0, H, 
B E3) EH E2]- This interest is because of its relevance 
to quantum gravity theory E3, EH and the success- 
ful fabrication of magnetic nanostructures with curved 
geometries QJ, El El El, E9, HI- In general, the fi- 
nite curvature of the underlying geometry may alter the 
geometric symmetries of the Hamiltonian describing the 
system; this alteration results in a qualitative change in 
the critical properties of the system. In fact, it has been 
reported [H, [H, [Hj] that such curvature-induced alter- 
ations occur when the Ising lattice is assi gne d to a surface 
with a constant negative curvature [25l l26j|. Significant 
shifts in static and dynamic critical exponents toward the 
mean- field values were unveiled in Refs. [HI, HH , and the 
mean-field property of the system was analytically proven 
by employing the renormalization group method [2J] . 

The abovementioncd findings led to the exploration of 
surface curvature effects in various kinds of statistical lat- 
tice models m m, m m m m, m ii m, m H3, m, 

HH, USE)]. For instance, the g-state Potts model on the 
negatively curved surface exhibited a first-order phase 
transition when q > 3 35], and the XY model on the 
same surface showed the absence of Kosterlitz-Thouless 
transition due to strong spin- wave fluctuation [3lT |. It 
is noteworthy that the concept of negatively curved sur- 
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faces (or spaces) is also relevant in many fields where the 
geometric character underlying the system is ofgreat sig- 
nificance, such asglass science [H, EE EE EE EE EE 14a . 
alasma ph ysics lip. l50l. 5ll . |52| , quantum transport [53l . 

string theory [Bl^ , and cosmol- 



|54[, chaos 15 
ogy [H M 

Recently, Baek et al. [39( studied the percolation tran- 
sition in negatively curved lattices. They found two dis- 
tinct percolation thresholds - one that corresponds to the 
occurrence of a single infinite-sized cluster and the other 
that corresponds to the occurrence of many infinite-sized 
clusters connecting a site deeply inside the lattice to that 
lying at the outmost boundary of the lattice. The latter 
threshold originates from an exponential increase in the 
total number of sites, N cx e L , with the linear dimension 
of the lattice, L (for the definition of L, see Section II in 
the present article). The exponential increase in N leads 
to an important consequence: the ratio of the perimeter 
to the area of the lattice remains finite even within the 
thermodynamic limit (L — > oo). This result is in con- 
trast with the case to that in the case of a flat plane, 
wherein N cx L 2 , i.e., the ratio becomes zero for large 
values of L. Similar percolation thresholds have been ob- 
served in enhanced binary trees [62j whose lattice struc- 
ture is quite analogous to that considered in Ref . [3§] . It 
is conjectured that the two-stage percolation transitions 
indicate the non-uniform ordering of the corresponding 
Ising model; that is, non- vanishing boundary effects even 
for large ./V may cause spatially non-uniform growth of 
ferromagnetic Ising domains under cooling, wherein the 
ordering process near the outmost boundary differs from 
that deep within the lattice. The precursor of such non- 
uniform growth was discussed in Ref. [9( ; the stable Grif- 
fiths phase was found near the outmost boundary at a 
temperature lower than the mean-field transition tem- 
perature. 

In the present article, we examine the ordering pro- 
cess of the negatively curved Ising model at temperatures 
close to and far greater than the mean-field transition 
temperature in order to clarify the boundary-spin contri- 
butions to the domain growth and domain distribution 
during cooling. Using Monte Carlo simulations, we have 
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FIG. 1: Regular heptagonal lattice established on the 
Poincare disk. Here, the number of concentric layers of hep- 
tagons L is 3, and the total number of sites is 112. All hep- 
tagons depicted within the circle are congruent with respect 
to the metric given in Eq. Im- 



proved that short-range order exists at temperatures far 
greater than the transition temperature. The existence 
of short-range order at such high temperatures is pecu- 
liar to the Ising model of negatively curved surfaces and 
is thus an important consequence of non-zero boundary 
effects that arc observed even in large systems. 



II. REGULAR LATTICE ON A NEGATIVELY 
CURVED SURFACE 

This section gives a brief summary on the construction 
of regular lattices on a surface with negative curvature. A 
surface with constant negative curvature can be defined 
as a single sheet of a two-sheeted hyperboloid expressed 
by 



y 



-i (*>i), 



(i) 



which is constructed in the Minkowski space endowed 
with the Minkowskian metric ds 2 = dx 2 + dy 2 — dz 2 . 
Although this definition is exact, it is inconvenient for 
computations since three coordinates are used to describe 
a geometry that has only two degrees of freedom. We 
thus use an alternative representation of the surface - 
the Poincare disk representation - that is obtained by 
projecting the upper hyperboloid sheet onto the x-y plane 
by using the following mapping: 



(x,y,z) 



y 



l + z 1 + z 



(2) 



As a result of this mapping, the upper hyperbolic sheet is 
transformed into a unit circle on the x-y plane endowed 
with the metric 



ds 2 



f(dx 2 + dy 2 ), f 



(1 — x 2 — y 



2\2 ' 



(3) 



This unit circle is referred to as a Poincare disk and serves 
as a convenient representation of the surface with a con- 
stant negative curvature. The Gaussian curvature n on 




FIG. 2: Graphical representation of the ratio of the number 
of outmost boundary sites to the number of total sites, as 
expressed by B(L) = [N{L)-N(L-\))/N(L). The horizontal 
axis represents the number of concentric layers of heptagons 
L. Inset: Semi-logarithmic plot of N(L) vs. L. 



the disk, which is obtained by [23J 

__u&_ 

f \ dx 2 dy 2 



log/, 



(4) 



takes the value of k = —1 at arbitrary points on the 
disk. The boundary of the disk corresponds to points at 
infinity in the hyperbolic plane. 

An infinite variety of regular polygonal lattices can be 
built on the Poincare disk. All the lattices satisfy the re- 
lation (p — 2)(q — 2) > 4, where p is the number of edges 
of each polygon, and q is the number of polygons around 
each vertex [26|. The heptagonal lattice of {p, q} = {7, 3} 
is considered in the present work. Figure[T]shows the hep- 
tagonal lattice represented as a Poincare disk. Although 
polygons depicted in the figure appear to be distorted, 
all of them are exactly congruent with the metric given 
in Eq. ([3]). The size of the entire lattice is characterized 
by the number of concentric layers of heptagons, denoted 
by L, which effectively serves as the linear dimension of 
the lattice. 

For a given L, the total number of sites AT is expressed 
by 



N(L) 



7 for L = 1, 

7 + 7J2 [u+ (v+f + u- («_)' 



for L > 2. 



(5) 



where u± — 2 ± y/S and v± = (3 ± v5)/2; refer to Ap- 
pendix for the derivation of Eq. |[5|). The inset of Fig. [2] 
shows the semi-logarithmic plot of N(L) with L. The 
plot shows that N increases exponentially for L ^> 1. In 
fact, we obtain 



B(L) 



N(L) — N(L — 1) 



1-— , (£->oo) (6) 



FIG. 3: (Color online) (a) Double logarithmic plot of spontaneous magnetization m(T) (thin line), and two-spin correlation 
function C{rij) (thick lines), (b) Correlation functions of two spins both of which lie within an Lth concentric layer. 



which implies the boundary contribution quantified by 
-B(-Zj) does not become zero, but remains finite within 
the limit L — > oo. It is emphasized that non- vanishing 
property of B(L) as well as the exponential increase in 
N(L) with L is a result of the constant negative curvature 
of the underlying geometry, and it is the reason behind 
the non-trivial critical behavior of embedded lattices, as 
proved in previous studies (22l . [231 ] . 



III. NUMERICAL METHOD 

The current study aims at exploring the ordering mech- 
anism of the hcptagonal Ising model with ferromagnetic 
interaction. The Hamiltonian of the heptagonal Ising 
model is given by 

H = -J^s iSj , 8i = ±l, (7) 

where denotes a pair of nearest-neighbor sites on 

a lattice with free boundary conditions. The coupling 
strength J is a constant, considering the fact that all 
the Ising spins are equally-spaced on the Poincare disk. 
Throughout this paper, J/ks and J are used as the units 
of temperature and energy, respectively. We have em- 
ployed Monte Carlo simulations to calculate the order 
parameter m = (sj) and the two-point correlation func- 
tion C(rij) = (si ■ Sj) with = |r$ — rj\. Here, 
represents the location of the ith spin on the Poincare 
disk, and the distance ry is a measure of the number of 
bonds along the shortest path between the ith and jth 
sites. Configuration space sampling has been carried out 
by using a cluster-flip algorithm [631 ] . and averages have 
been calculated by considering 10 5 samples. 



IV. RESULTS AND DISCUSSION 

Figure [3] shows the temperature dependence of 
C(rij,T) for different r^'s. We have set L — 9 (i.e., 
TV = 40432) in all calculations. We found that for r v] < 4, 
C(rij) shows power- law decay with T and thus survive 
at T > T c ~ 1.25 [12] though m(T) is almost zero there. 
These survival of Cfaj) indicates that small-sized ferro- 
magnetic domains remain active even at high tempera- 
tures, although they do not contribute to m(T) because 
of the change in sign between positive (sj = +1) and 
negative (si = —1) domains. To explore the variety 
of domain sizes and examine boundary effects on size 
variations, we evaluate the correlation function, which is 
denoted by Cl (rtj , T) , of two spins that lie in the Lth 
layer. Figure G^b) shows the T-dependence of Cl for 
Tij = 1,2,4, in which L is varied from 6 to 9 (data for 
L < 5 are omitted since they are indistinguishable from 
data for L = 6). Interestingly, the data of Cl for each 
Tij value collapses onto a single curve at T > T c . This 
collapse implies that the domains that continue to exist 
at T > T c are of the same size and are uniformly scat- 
tered over the entire lattice, regardless of their distance 
from the outmost boundary. 

It is emphasized that the persistence of domains at 
such high temperatures is supplementary to previous 
findings [H, [22j related to low-temperature behaviors of 
an identical system. It was found that at T < T c , the 
central region of the lattice (far from the outmost bound- 
ary) exhibits an ordered phase that is a consequence of a 
mean- field transition [22j , while the outer region (close to 
the boundary) exhibits a stabilized Griffiths phase that is 
free from extrinsic disorder [§]. Note that the former re- 
sult indicates the presence of a paramagnetic phase in the 
central region at T > T c , where Ising spins are randomly 
oriented due to large thermal fluctuations. Hence, one 
would expect the formation of ferromagnetic domains to 
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FIG. 4: (Color online) Spontaneous magnetization mi(T) 
within a circular region T>l enclosed by the Lth layer. 

be hindered at T > T c ; however, contrary to the expec- 
tation, the formation is observed in the current simula- 
tions. This apparent contradiction is attributed to the 
finite-size effect. For a moderagely large L, the strong 
boundary-spin contributions penetrate to the center of 
the lattice, as a result of which small-sized domains arise 
not only near the boundary but also in the central region. 
We conjecture that the domain in the central region dis- 
appears if a sufficiently large value of L is considered: 
however, huge computational costs are involved when the 
value of L is large enough. 

To gain a better insight into the above issue, we eval- 
uate the quantity i71l{T) = jjJ2 s ( r i c where T>l 
is a circular region enclosed by the Lth layer. Figure 2] 
shows the T-dependence of (T) for different values of 
L. It is found that above T c , m/,(T) for L = 3 is finite 
whereas m(T) becomes almost zero. This indicates an 
incomplete cancellation between the positive and nega- 
tive domains in the region Dl=3 due to its small size. 
We emphasize that ra^T) should become zero if T>l is 
in the paramagnetic phase with no finite-sized domains. 
Therefore, the finiteness of rriL=-i{T) at T > T c is an- 
other piece of evidence that supports the persistence of 
small-sized domains. We have also confirmed that finite 
rriL=3(T) at high T tends to disappear with an increase in 
L, since boundary-spin contributions are prevented from 
penetrating to the center of the lattice. 

Figure [5] illustrates the ordering process of our mod- 
erately large Ising lattice with a decrease in the temper- 
ature. At T T c , small-sized ferromagnetic domains 
(indicated by open and solid circles) are randomly em- 
bedded in the paramagnetic phase which is shown in gray. 
These domains are distributed homogeneously across the 
lattice, but give rise to finite mt(T) for small values of L 
due to the incomplete cancellation between the positive 
and negative domains. Subsequently, a decrease in T re- 
sults in the growth of domains within the inner region. 
Eventually, the ferromagnetic phase is formed through 




FIG. 5: Illustration of the ordering process: the left panel 
represents T < T c ~ 1.25; the middle panel, T ~ T c ; and the 
right panel, T > T c . The open and solid circles indicate fer- 
romagnetic domains with positive and negative directions, re- 
spectively. The gray area represents the paramagnetic phase. 

the mean-field phase transition. Nevertheless, near the 
outmost boundary, small domains remain active and fluc- 
tuate on a large time scale, as observed in Ref. @. 

V. CONCLUSION 

In the present work, we have considered the ordering 
mechanism of the Ising lattice model assigned to a neg- 
atively curved surface. We have found that small-sized 
domains survive at temperatures much higher than T c , 
at which the inner bulk region of the lattice goes through 
mean-field phase transition. The existence of small do- 
mains is attributed to the non-zero boundary-spin con- 
tribution that is unique to negatively curved surfaces. 
Our results are consistent with those of previous studies: 
the Griffiths phase is present near the outmost bound- 
ary of the lattice, while the inner core region undergoes 
the mean-field phase transition. This consistency sheds 
light on the nature of phase ordering in a wide variety 
of statistical lattice models assigned to negatively curved 
surfaces, most of which need to be further studied. 
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APPENDIX A: PROOF OF THE FORMULA (5) 

This Appendix is devoted to the derivation of Eq. ((SJ) 
- the explicit form of N(L) as a function of L. We 
classify all the congruent heptagons in the lattice into 
two groups: (i) heptagons having a common edge with 
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the adjacent heptagon lying in the inner concentric layer 
(referred to as "A" -heptagons) and (ii) those having two 
common edges with the adjacent one heptagon in the 
inner layer (referred to as "B" -heptagons). Figure HJa) 
illustrates our classification, in which either of the sym- 
bols A or B is assigned to all heptagons except for the 
central one. 

We see in Fig. [5] that the second innermost layer 
(L = 2) consists of seven B's, the third innermost layer 
(L = 3) consists of 14 B's and seven A's, and so on. We 
identify the four sites within a hatched region as those 
attached to the i3-heptagon; in a similar way, the three 
sites enclosed by a curve are identified as those attached 
to the ^4-heptagon. Then, we can say that the four sites 
of a i?-heptagon located at L — 2 "engenders" one A- 
heptagon and B-heptagons at L = 3. Similarly, it follows 
that an A-heptagon at L = 3 generates a pair containing 
an ^4-heptagon and a L?-heptagon at L = 4 (although it 
is not shown in Fig. [B^a)). This proliferation process 
is summarized by the diagram in Fig. [5Ja); each A (or 
B) block symbolizes an ensemble of three (or four) sites, 
and a block in the Lth layer creates several blocks in the 
(L + 1) layer. The diagram illustrates a method to count 



the total number of sites contained in a given Lth layer. 

We now have all the ingredients to derive the formula. 
Let us denote the number of A- (or B-) heptagons in the 
Lth layer by (or bi) and set et2 = 0, 62 = 7. It follows 
from the diagram that for any L > 3, 



a>h = a-L-i + b L -i, b L = a L - 1 + 2b L - 1 . 



(Al) 



We eliminate &l from Eq. (|Aip to obtain — 3a£_i + 
a.L-2 = whose solutions under the conditions a 2 = 
and 03 = 7 are given by 



7V5 



L-2 



,L-2 



for L > 3 (A2) 



where v± = (3± V5)/2. Hence, we have for L > 2, 
7a/5 



br = 



10 



L-2 



(A3) 



where w± = ±1 + \fb. As a result, the desired formula 
Eq. (0 - for N(L) = 3a L + 4b L can be obtained. 
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